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Abstract 

We study production of sterile neutrinos in the atmosphere and their detection at 
Super-Kamiokande. A sterile neutrino in the mass range 1 MeV < Mjv ^ 105 MeV 
is produced by muon or pion decay, and decays to an electron-positron pair and 
an active neutrino. Such a decay of the sterile neutrino leaves two electron-like 
Cherenkov rings in the detector. We estimate the sterile neutrino flux from the 
well-established active neutrino fluxes and study the number of the decay events 
in the detector. The upper bounds for the active-sterile mixings are obtained by 
comparing the 2e-like events from the sterile neutrino decays and the observed 
data by Super-Kamiokande. The upper bound for the muon type mixing 0^ is 
found to be |G/,|2 < 5 X lO'^ for 20 MeV < Mat < 80 MeV, which is significantly 
loosened compared to the previous estimation. We demonstrate that the opening 
angle and the total energy of the rings may serve as diagnostic tools to discover the 
sterile neutrinos in further data accumulation and future upgraded facilities. The 
directional asymmetry of the events is a sensitive measure of the diminishment of 
the sterile neutrino flux due to the decays on the way to the detector. 



1 Introduction 



The existence of nonzero neutrino masses has been confirmed in the last few decades 
and stimulated the activities aiming for the theory beyond the Standard Model. Among 
various possible ways to introduce the neutrino masses into the model, adding gauge 
singlet fermions (right-handed neutrinos) is one of the most economical and attractive 
methods. In particular, the seesaw mechanism [1] naturally accounts for the observation 
that the neutrino masses are very small compared to the other fermions. Moreover, it also 
gives a natural prescription for the baryon asymmetry of the universe by the leptogenesis 
scenario |2]. 

Besides the usual higher mass scales of the right-handed neutrinos around the grand 
unification scale, the lower mass ranges are also interesting and rich in phenomenology. 
For example, a keV sterile neutrino is a viable dark matter candidate [21 H] and accounts for 
the pulsar velocities [5] due to feebleness of its interactions. Two quasi-degenerate sterile 
neutrinos in the mass range (9(10^^) — (9(10) GeV provide viable baryogenesis scenario 
alternative to the leptogenesis [HI [TJ [H] . Remarkably, these excellent features originate in 
a single framework so called i/MSM [H [7] , which is an extension of the Standard Model 
with three generations of the right-handed neutrinos. Due to lower threshold energies of 
the production, such "light" sterile neutrinos are more likely to be tested in existing and 
forthcoming experiments than the usual right-handed neutrinos with super heavy masses 
(for example, see Ref. and references therein). Mass range by mass range, the sterile 
neutrinos may provide unique signals in various circumstances. It is thus interesting to 
study how to produce and detect the sterile neutrinos. 

In this paper, we focus on the sterile neutrinos with the mass range 1 MeV < Mjv ^ 
105 MeV and study their production in the atmosphere. The previous neutrino exper- 
iments, including peak searches in the meson decays [TOl [TTl [121 UHl IHj and the decay 
search with accelerators [l5l[T6], have placed certain bounds on the active-sterile mixing 
matrix in this mass range. The atmospheric sterile neutrino provides an independent and 
complemental test to these experiments with artificial neutrino sources. In particular, the 
peak searches of the pion and kaon decays put no stringent bound on the muon type mix- 
ing angle around = 40 MeV. This spot will be effectively probed by the atmospheric 
sterile neutrino since the sterile neutrinos in this mass range can be copiously produced 
by the muon decays similar to the atmospheric active neutrinos. 

The production of sterile neutrinos in the atmosphere and their detection at Super- 
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Kamiokande (SK) has been discussed in Ref. [T7]. In this work, we extend and improve 
their work in both sterile neutrino production and its detection. In the flux estimation, we 
carefully consider the production of the sterile neutrino not only by the decay, but also 
by the tt^ decay. The sterile neutrino flux does not receive the vr^ contributions for ^ 
35 MeV when the mixing with electron is sufficiently small (as we will show later, this is 
indeed the case by considering the bounds on the mixings from direct search experiments) , 
being different from the and z/^ fluxes. We also take into account the energy distribution 
of the daughter (the sterile neutrino) in the parent (mainly muon) decay and the phase 
space suppression due to the sterile neutrino mass. With the sterile neutrino mass around 
the neighborhood of the muon mass threshold, the flux gets suppressed and the mixing 
bounds are signiflcantly loosened compared to Ref. [IT]. Furthermore, we study the 
kinematics of the detection process N ^ e~e^v m detail, which include the distributions 
of the opening angle and the visible energy of the emitted e~ and e"*". We estimate the 
upper bounds on the mixing angles by using 1489 days SK data [TS^. In this estimation, we 
shall apply the cut to the events by requiring the visible energy being larger than 30 MeV 
in accordance with the SK data. This treatment also changes the previous results. 

This paper is organized as follows. In Section [2], we calculate the sterile neutrino fluxes 
produced in the atmosphere. In Section [3], the sterile neutrino decay — )■ e'e^v and its 
event rate at SK are studied. Section H] is devoted to conclusions. 

2 Fluxes of the sterile neutrinos 

We consider a gauge-singlet fermion with the mass range 1 MeV < M^v ^ 105 MeV 
which mixes with the left-handed neutrino {a = e, fi, r) as 

Z/q = f/m Vi + QaN, (1) 

where Uai is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, Vi {i = 1,2,3) are 
the mass eigenstates of the active neutrinos. The parameter 0^ is the mixing between 
active and sterile neutrinos, which rules the interaction strength of A^. The extension to 
the multi-generation case is trivially done by replacing 6q,A^ with QaiNj. Throughout 
this work, we focus on the case where the sterile neutrino is mainly mixed with z/^ for 
the shake of simphcity. That is, we assume |9^| ^ |9e|, |0t|, unless otherwise stated. 

Indeed, the electron type mixing Be is much more severely constrained than the other 
two parameters. The peak search with vr^ — )■ e^A^ mode suggests |6eP ^ 10~^ — 10~^ for 
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40 MeV < Mtv ^ 140 MeV, while the decay search with accelerators indicate the upper 
bound for |e^p varies from 10"^ to 10"^ in the mass rage 10 MeV < M^v < 100 MeV [H]. 
In the case of dominance, the main decay mode of is — > 3z/ and the subdominant 
mode is N ^ e^e^u^ conducted by the neutral currents. The decay width of each process 
is given by 

r(iv ^ 3.) = ^£^1^, (2) 

r(A^ ^ e~e+z/) = r(A^ ^ 3u) Q - sin^ + 2 sin^ 9w^ ■ (3) 

In this paper, we focus on the case where the neutrinos are the Majorana particles. 
Then the lifetime of the sterile neutrino is given by r ~ l/2r{N — Su) = 1.1 x 
10-^^ (^^li^) (s) and the corresponding decay length is ~ 0.33 ^ (^^i^) (km) 
for = 100 MeV sterile neutrinos. For = 10~^ and = 100 MeV for example, 

the decay length is about half of the earth radius ~ 6,400 km. The branching ratio 
of the detectable mode is ~ 1/4 — sin^ 9]v + 2 sin^ 9w = 0.13. 

The sterile neutrinos are produced from charged pions and muons in the same manner 
as the active neutrinos. The main production processes are 

7r± ^ fi^ N, 

The former channel is open only for Mn < m.,^± — ^ 35 MeV. The kaon contribution 
is negligible in the energy range of present concern. In the muon neutrino production 
for example, K"^ contribution becomes significant only in the high-energy regime > 
100 GeV [ig. 

In fact, the vital part of the sterile neutrino flux is low-gamma regime 7 ^ 10 for 
an evaluation of the detectability of the atmospheric sterile neutrinoi*]- Namely, we are 
interested in the sterile neutrino spectrum at most up to 1 GeV for the mass range of 
1 MeV < Mtv < 105 MeV. As we will discuss in Section |3] in detail, the sterile neutrinos 
are detected by the decay A^ — > e^e'^u which leaves 2e-like rings in the detector. If the 
gamma factor of the sterile neutrino is too large, the two fuzzy rings will overlap each other 
and the separation of the rings becomes difficult. Indeed, it turns out that 50% (80%) of 
the whole events will be observed with the opening angle less then 20° for 7 = 6 (10). 



*Here and henceforth, we use the symbol 7 to denote the gamma factor of the sterile neutrino unless 
otherwise mentioned. 
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It is in general a complicated task to calculate the flux for such low-energy regime since 
it is affected by solar activity and the geomagnetic field. In this work, we do not compute 
the flux from the primary cosmic ray spectrum directly, but we try to reconstruct a 
reasonable parent's (muons and pions) spectrum from the well-established active neutrino 
fluxes available in literature [201 EI] and then evaluate the sterile neutrino flux from the 
reconstructed parent fluxes. 

The procedure goes as follows. We start from the and z/g fluxes at SK site ^20j to 
estimate the muon fluxes. The evolution of the neutrino flux (j)p^{E^, t) is described by [19] 



dME.,t) _ r,E^^'-^iE,E.l (4) 



dt Je^ X{E,t) TdE, 

where t is the slant depth, (j)^-{E,t) is the muon flux, E^ (E) is the neutrino (parent fi~) 
energy, A = p{t)^{E/m^y - 1/r is the muon decay length multiplied by the density of 
the atmosphere p{t), T is the decay width for fi^ — )■ e~z/eZ/^ in the laboratory frame. By 
integrating (j4]) with respect to t, one finds 

ME.,tf) = JjE <D,-(E) f^^iE,E^), V(^) = jy^^^^, (5) 

where tf is the depth corresponding to the sea level. The sterile neutrino flux (f)fyf{EN,t) 
obeys the similar equation to the active neutrinos if the sterile neutrino decay in the 
atmosphere is negligible. In fact, the decay length for — )■ 3z/ is much larger than 
the altitude of the mesosphere edge ~ 100 km for most of the parameter space. Only 
exception is the regime |9^p > 6.6 x 10"^ ( ^"a/^''^ ) ( looMev ) where N Su is so rapid 
that some of the sterile neutrinos do not reach the detector. This means for example, the 
following estimation of the sterile neutrino spectrum is not reliable and may get additional 
suppression for E^ < 100 MeV in the parameter regime IQ^]"^ > 6.6 x lO^'^ {j^l^^^ ■ 
Putting aside this strong interacting regime, we have 

where A' and V are the decay length and width for the mode p~ — )■ e~UeN. By using 
A' = (r/r')A, one finds 



{E^,tf) = r <I>,-{E)-^{E,E^), (7) 



En 

T' M 
— = |e,P ( 1 - + 8r^ - - 24r^ln(r) ) , r = — 
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Here the integrated muon flux [E) is obtained by fitting the left-hand side of ([5]) 
with an assumption of power-low behavior of ^^-{E). The same exercise is applied for 
i'e-fJ''^-N chain and then the flux for the Majorana is obtained as the sum of these 
two contributions. Details for the energy distributions in muon decay are presented in 
Appendix lA.li 

In the mass regime < m^± —m^, pions also contribute to the sterile neutrino flux. 
Roughly speaking, vr^ contribution is the same as the muon contribution, so that the flux 
of the sterile neutrino in M^r < wi,r± ~ regime is twice as large as the higher mass 
regime up to the threshold effect. The vr^ contribution is estimated in the same way as 
the muon. For the vr^ case, however, we need z/^, z/^ fluxes of vr^ origin to reconstruct the 
integrated vr^ fluxes from the muon neutrino input. We assess the z/^, of tt"^ origin by 
subtracting the z/^, z/^ of /i^ origin from the full z/^, z/^ data in Ref. |20j, where t'/^, z/^ of 
^± origin is calculated from the integrated muon fluxes <l>^±(£') reconstructed via Eq. (I5l). 
Finally, the sterile neutrino flux of tt^ origin is calculated by 

iEr,,tf) = / dE^ME)-—^{E,E^), (9) 

where F^^ and F^ are the decay width for vr^ — )■ /i^z/^(z/^) and n"^ — )■ fi'^N, respectively. 
Here the branching ratio is given by 

T^f I (1^2 _|_ iy,2i { r'^ T*^^^ 

- le W l-2(r^ +r^) + (r^ - r^)^ ^1^ >L 

F^ ~ ' V ^VN + r^l) + yrN ^^i) ^2(^1 _^2y 

Tn = , r^, = — ^. (10) 

Fig. [1] shows the sterile neutrino fluxes for several examples of M^r. The six types of the 
lines are for = 3, 10, 30, 40, 60, 80MeV, respectively (from the left). It is seen that 
a naive estimation (pN = 1©^^^'^'^^+!^^ made in Ref. [Hj is good for the energies above the 
mass threshold with Mjv < 30 MeV, but it breaks down significantly for Mn > 30 MeV 
in such a way that the calculated fluxes are suppressed as Mn increases. This rapid fall 
off with respect to is due to the phase space suppression of the muon decay (See 
Eq. ([8])). On the other hand, the phase space suppression of n"^ decay is significant only 
at the vicinity of the threshold m^± —m^. As increases, the flux thus suddenly drops 
by factor of two at the threshold ~ 35 MeV above which the tt"^ contribution vanishes. 
The sterile neutrino fluxes presented here are averaged over the zenith angle and given at 
the middle level of solar activity following the input neutrino fluxes in Ref. [20] . 

The uncertainty of the reconstructed parent fluxes and $,r± is small for high- 
energy regime, but it becomes larger as the energy is lowered so that the resultant sterile 
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Figure 1: Fluxes of the sterile neutrinos. The six curves are for = 
3, 10, 30, 40, 60, 80MeV with the mixing = 10~^. The triangles and squares show 
the active neutrino fluxes from Ref. [20j. The black circles show |B^p0;^^+p^^. 

neutrino fluxes carry certain ambiguity around En — M^. This is because high-energy 
part of the active neutrino flux is well fitted by a single-power low, while low-energy 
part are complex and many choices are available as fitting function. In this analysis, we 
have used the fitting function ^h^l/{^h + ^l), where and $x are the high and 
low-energy part of the fitting function. For ^h, we have taken = aul^" with some 
constants an-, bn and the gamma factor 7 of the parent particle of interest. For we 
have examined two options $l = aiP^^ or $i = a^'y^'^ , where (3 is the beta factor of 
the parent particle. While with the former the flux rises from the mass threshold 
and continuously shifts to constant behavior and is finally reduced to ^h, the latter $2. 
has discontinuity at the mass threshold and the peak is slightly shifted to lower energies 
compared with the former one. It turns out that two different fitting bring at most 10% 
difference in the total number of events. We would like to emphasize that the salient 
feature of the phase space suppression is however independent from the fitting scheme 
and the conclusions derived from this effect are valid in what follows. The fluxes shown 
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in Fig. [Dare calculated with the option $l = aL(3^^ and these fluxes are used in the 
following analyses. 

3 Decay of the sterile neutrinos and the event rates 

The sterile neutrinos are produced in the atmosphere and reach the SK site to leave 
2e-like events via — )■ e~e"'"z/ and — > e~^e~u. Suppose that the flducial volume of SK 
is represented by a sphere of radius r. The event rate is then given by 

roo r i'l+2r \ i' f 1 //F 

Rate = dEN r^dn(t)N / dl'—e~~d dX—-^{EN,X), (11) 
Jmm J Ji J Tat dX 

where and are the decay width and length for the signal decay of A^, X is an 
observable of interest, e.g., the invariant mass of the momentums of e~ and e"*", and 
/ is the flight distance of the sterile neutrinos. The total decay length is given by 
Ad = AdA^/(Ad + A'^) where A'^ is the decay length for A^ — )■ 3z/. Since r is much smaller 
than the earth radius i?®, the flight distance / is well approximated by 

/ = / ° for < < I, 

\ -2i?eCose for f < ^ < vr, ^ ' 

in terms of the zenith angle 9. The total events are given by the integration over the solid 
angle J dVt = J^^ dcosO J^^ d(f). By integrating over all possible flnal states of e~e~^iy and 
taking A^ ^ 2r, Eq. ( fTTj) becomes 



Rate = 3Vm / dEj^ -7-<Pn 

Mm Ad 



2Ra 



2Rq 



+ 1 



(13) 



where Vfid is the flducial volume of the SK tank 22.5 kton which corresponds to 2.25 x 
10^ m^. Here the flrst (second) term in the square brackets is the contribution from the 
up (down)-going events. For A^^ ^ 2_R^, one flnds 

Rate = 6Vm / dE^ t-^tv, (14) 

JMm 

which is independent from the total decay width and depends only on the partial decay 
width Ftv responsible for the detection channel. The events take place isotropically due to 
the condition A^ ^ 2i?® with which the flight distance / becomes irrelevant to the event 
number. While if A^ < 2i?0 so that the sterile neutrinos partially vanish on the way to the 
detector, downward-going events dominate over upward-going event. The measurements 
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Figure 2: Left panel: Number of events (A^ — )■ e~e~^i') at SK per year. The red solid (blue 
dotted) lines shows the case when IG^^P = 10^^(10^^). The thick and thin lines are with 
and without the effect of decay from atmosphere to detector. Right panel: Up-down 
asymmetry in number of events. Here we take Mn = 60 MeV. 



of the direction of the total momentum of e~ and may thus clarify the indication of 
iV — )■ 3z/ decay taking place on the way to the detector. 

The left panel of Fig. [2] shows the total number of events per year as a function of 
the sterile neutrino mass M^. The red solid (blue dotted) lines shows the case when 
|0^P = 10~^(10~'*). It is seen that the event number is proportional to |0^|^. The thick 
and thin lines are with and without the effect of N decay from atmosphere to the detector. 
The thick lines are suppressed compared to the thin lines by the amount that the up-going 
events are reduced by the decay effect. The decay effect is negligible for |6^p = 10~^ 
while it makes certain difference for |0^P = 10~^. The right panel shows the ratio of the 
up-going events to the down-going events. As we will see later, the asymmetry is unlikely 
to be observed for |0^P ^ 3 x 10~^ and = 60 MeV where the number of events does 
not exceed the 2e-like ring data of SK. 

If Gr is switched on, however, it hastens both decay N ^ Su and A^ — )■ e'e'^u while 
the production processes are kept unchanged. With being finite, the decay width of 
each mode is obtained by replacing |6^p with |6^p + |6rP in (12]) and ([3]). By virtue 
of 0e = 0, the above two decay modes are induced only by the neutral current so that 
the effect of 0,- 7^ simply appears as such a simple replacement. On the other hand, 
the production processes of A^ are not affected by 0^ since the decays of vr^ and /i^ are 
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induced by the charged current and cannot involve a tauon. Under the assumption that 
9e = 0, the production processes thus involve only 0^. In the flux equation IQ, the decay 
term of would not be negligible if |Gt-| is large. This effect however gives negligible 
contribution to the ratio of upward to downward-going events since the thickness of the 
atmosphere is small compared to the earth's diameter. 

An interesting possibility in view of the up-down asymmetry is therefore |6t-| ^ |0^| 
with which the flux is feeble while the decay probability is high. In such case the event 
rate is proportional to |0^p| 10,-^ instead of |0^|^, where comes from the production 
and |0rP is from the detection. Let us set |0^p||0T-p ~ 10~^ to obtain C(10) events 
per year. Then, for instance, |0tP ~ 10~^ and |0^P ~ 10"'' lead to a clear asymmetry 
jyup/jydown _ q^^^ (Notlcc that the right panel of Fig. [2] can be read plot for 1 0^1 in 
the case where |0r| ^ asymmetry is observed, it may indicate a hierarchical 

structure of \Qr \ ^ 

The master formula for the event rate (fTTl) involves the product of the production and 
the detection probabilities. That is, the product of the squared modulus of the production 
and the detection amplitudes. This means that (fTTj) is valid in the case where the coher- 
ence between the propagating neutrino states is lost. Since we have the sterile neutrino 
and the active neutrinos in the theory, one may wonder if the oscillation between the ac- 
tive and sterile states occurs. However, the coherence between the active and sterile states 
is explicitly violated at the detection point; the sterile neutrino decays while the active 
neutrinos do not. Furthermore, even if the coherence were not violated by the detection, 
such an oscillation would be so rapid that the effect is averaged out, which also insures 
the validity of using (ITT!) . In fact, with = lOMeV and a typical energy of 100 MeV 
for example, the oscillation length L""^^ would be L""^^ ~ ( ioomcv ) ( ^° Am°^^ ) 10"^^ 
the present setup, it is quite difficult to observe oscillatory behavior between the active 
and the sterile neutrinos. 

3.1 Invariant mass and bounds for the mixing 

The sterile neutrino decay — )■ e'e'^u in the SK detector produces two fuzzy Cherenkov 
rings. A possible background for this signal is tt" — 27, where tc^ are mainly created via 
the neutral current interactions of the atmospheric active neutrinos. The gamma rays 
in the final state develop into electromagnetic showers and create two fuzzy rings which 
cannot be distinguished from the signals. 

The above two process, however, leave different signatures in the Lorentz invariant 
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Figure 3: Invariant mass distribution of the electron pairs from the sterile neutrino decay 
— )■ e~e+z/. The solid (red), dashed (green), and short-dashed (blue) lines are for 
Mm = 80, 60, 40 MeV with the mixing |0^P = 10"'^. The points with error bars are the 
data of fully-contained 2e-like rings from Ref. [T8] . 

mass (squared) of the 2e-like rings. For the sterile neutrino event — )■ e~e^z/, it is given 
by Mgg = {pi +P2Y, where Pi{p2) are the four-momenta of e~{e~^). In the rest frame of 
the sterile neutrino, it reads 

Ml ^ M%- 2MnE,, (15) 

where E^, is the invisible energy carried away by the active neutrino. Since the neutrino 
carries about 1/3 of the parent energy on average, M^ from the sterile neutrino events 
follows (Mgg) ~ Mff/3, while the invariant mass squared of the two photons by vr^ decay 
sharply peaks at ttl'^o. By measuring the invariant mass distribution, we may obtain 
information on the sterile neutrino mass. 

It is more practical to work with M^e than M^. The Mge distribution is given by (See 
Appendix IA.3I for details) 



where Uee = Mee/M^- The distribution ( fT6l) has the maximum at Uee ~ 0.5, and the 
averaged invariant mass is given by {M^e) = J^ ^'^ dM^e^-^^ = 0.508 M^. Note that the 
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Figure 4: Visible energy distributions for M^r 
Here we take |6^P = 10~^. 



20 MeV (left) and Mn = 60 MeV (right). 



distribution function ( lT6l) also holds in the laboratory frame due to the Lorentz invariance 
of Mgg. We calculate the event rate for each invariant mass by using ( IT6l) in ( ITTll . 

Fig. [3] shows the number of events in each invariant mass bin. The solid (red), dashed 
(green), and short-dashed (blue) lines are for = 80, 60, 40 MeV with the mixing 
|0^P = 10~^. The points with error bars are the 1489 days data of fully-contained 2e-like 
rings from Ref. [18]. Following ( 1T6|) . the events are distributed to the range < Mee < 
and frequently seen at ~ 0.5Mn. From = 40 MeV to = 60 MeV, the hight of 
the peak grows although the possible range of M^e, to which the events are distributed, 
is broadened. This is because the total number of event increases from = 40 MeV 
to Mat = 60 MeV (see Fig. [2]). On the other hand, the hight of the peak rapidly falls off 
from Mn = 60 MeV to Mn = 80 MeV since not only the range of Mee is broadened but 
also the total number of event decreases. 

A naive upper bound for 0^ would be determined by comparing the hight of each bin 
with the corresponding background for fixed values of Mn- However, it should be noticed 
that the SK data shown in Fig. [3] is for visible energy greater than 30 MeV, where the 
visible energy means the total energy of the two rings (hereafter we denote it by Eee)- 
Hence we must cut the events whose Eee is less than 30 MeV to make such a comparison. 
In addition. In addition, the opening angle between two momenta of e~ and must be 
sufficiently large in order for the event to be identified as multi-ring events. Thus one must 
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Figure 5: Opening angle 9ee distributions for Mjv = 20MeV (left) and = 60MeV 
(right). Here we take \Qfj,\'^ = 10~^. 

study how E^e and the angle cut reduce the number of events for a reasonable estimation 
of the allowed parameter range. 

Roughly speaking, the two cut define an effective range of the sterile neutrino energy 
E^- for a given value of M^r. In particular, the visible-energy cut defines a lower limits 
of Ejy. In the three-body decay — )■ e~e~^i>, the electron and positron carry on average 
about 1/3 of the parent energy each, so that the events with E'^r ^ 45MeV more or 
less pass the cut. Accordingly, the number of event receives significant reduction for 
Mm < 45 MeV while it does not for Mjv ^ 45 MeV. Fig. H] shows two examples of Eee 
distributions. It is seen that about 60% of the whole event is dropped for Mjv = 20 MeV 
but the cut is insignificant for Mn = 60 MeV. 

On the other hand, the angle cut does not dramatically change the number of events. 
In Fig. we show the opening-angle distributions for Mjv = 20 MeV and M^- = 60 MeV, 
assuming angle resolution of 10°. It is clear that the angle cut does not have as much 
impact as the visible energy cut for both cases if the two rings with 6ee ^ 10°, or more 
conservatively, 9ee ^ 20°, can be separately identified at SK. We found that efficiencies 
of the angle cut are 0.65 for = 10 MeV, 0.8 for Mjv = 20 MeV and 0.85 for Mjv = 
30 — 100 MeV with 6'ee > 20°. The sterile neutrinos with smaller mass receives more 
reduction than larger mass since light sterile neutrinos are more energetic. It is seen in 
Fig. [1] that sterile neutrinos are mostly populated in 7 ~ 4 — 5 for Mjv = 20 MeV while 
7 ~ 1 - 2 for Mtv = 60 MeV. 
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Figure 6: Correlation between E^^ and 6^^ for = 20 MeV (left) and = 60 MeV 
(right). The colors in each bin shows log^gl'^^^ber of events). jO/^P = lO""^ is taken. 

We present correlations between Ege and 6ee in Fig. |6] for completeness. For a lighter 
mass Mtv = 20 MeV, most sterile neutrinos are energetic 7 ~ 4 — 5 and the events are 
concentrated on lower-left region. For the heavier mass Mn = 60 MeV, the peak is shifted 
toward larger 6ee region since there are many sterile neutrinos with 7 ~ 1, and 6ee tends 
to be close to ~ 120° expected from the decay at rest. We would like to emphasize that 
these plots may be useful to separate the signal from the background if the tendency of 
the background is different from the signal. 

Having details for the visible energy and the opening angle distribution, let us come 
back to the invariant mass and estimate the upper bound of the mixing element G^. Fig. [7] 
shows the invariant mass distributions of events with a kinematical cut. As is already 
mentioned, lighter-mass case receives more significant reduction by the cuts. We would 
like to stress that, while the visible energy cut E^e > 30 MeV is indeed applied to the SK 
data, the angle cut is not the actual one which is applied to the data. In fact, the multi- 
ring identification process performed in the SK analysis |25] is much more complicated 
than the selection of opening angles and such a thorough analysis is beyond the scope 
of this paper. Since we do not know the minimal opening angle appropriate for the SK 
detector and could not find suitable references for this value, here and henceforth we 
assume it to be 20° just for example. We estimate the upper bound for 9^ by searching 
for the maximal value of with which none of the signal events of each bin exceeds the 
central value of the data [18] (shown by dots with error bars in the figures) for each fixed 
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Figure 7: Invariant mass distributions of envents at SK for M^v = 20 MeV (left) and 
Mtv = 60 MeV (right). The light gray region with red dotted line is without cuts, while 
the dark gray region with red solid line is with 9(.e > 20° and E^e > 30 MeV. 

value of Mm- 

Fig. |8] shows the upper bound in comparison with the previous result in Ref. p^. The 
bound is changed in three ways. First of all, in the large-mass region ^ 50 MeV, 
the bound is significantly relaxed due to the phase space suppression of the muon decay 
which has not been taken into account in the previous analysis. Second, in the small- 
mass region ^ 50 MeV, the event cuts (mainly the visible energy cut Ef.f. > 30 MeV) 
reduce the number of events so that the bound is pushed up from the curve without the 
cuts. Finally, careful handling of the tt"^ and /i^ contributions for the sterile neutrino 
production makes a small dip at Mn = nij^i — since the number of events are reduced 
for Mat > m7r± — according to the absence of the t^"^ contributions. 

The light-element abundances predicted by the Big Bang Nucleosynthesis (BBN) is 
kept unspoiled if the lifetime of is short enough such that the sterile neutrinos are 
cleared away before the onset of the BBN. According to Ref. [22], the successful BBN 
requires 



for 10 MeV < M^v < 140 MeV. Here this limit is valid for the Dirac type sterile neutrino 
which mixes with only one active flavor. For the Majorana case, this limit becomes two 
times weaker since the total decay width of the Majorana particle involves an extra factor 
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Figure 8: Upper bounds on the mixing element The red sohd and blue dashed 

lines are bounds with and without cuts O^e > 20° and E^e > 30 MeV, respectively. The 
black dotted line shows the bound in Ref. pTTj. 

of two. Most of the allowed region in Fig. [8] is in fact under tension with the limit (fT7j) . 
However, this tension is relaxed if is turned on. Namely, the lifetime is shortened by 
increasing |0t-|- The point is that the production of the sterile neutrino in the atmosphere 
is not affected by |0t|- Hence, by taking |0,-| ^ |0/^|, one can set a short lifetime enough 
to avoid the BBN limit while keeping the atmospheric sterile neutrinos consistent with the 
SK data. Notice that Fig. [8] can be read as a plot for 1 6^ 1 1 ©t- | if | ©r | ^ 10^1- Interestingly, 
the directional asymmetry (see Fig. [2] and related discussion) is likely to be observed in 
such a circumstance. 

If the visible energy cut En > 30 MeV is switched off, the sterile neutrino events in 
the first three bins increase as shown in Fig. [TJ while the background and its Mote Carlo 
simulation might change differently depending on the source of the background. Such an 
analysis is important especially for light sterile neutrinos with Mjv < 30 MeV but beyond 
the scope of this paper. 
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Figure 9: Upper bounds on the mixing element |6^P (the red sohd curve) in comparison 
with vrious experiments. The dotted (purple) curve is from the peak search in pion 
decay [TT] (labeled by 7r+ — )■ ^'^N), the dot-dashed (light-blue) line is from the peak 
search in pion decay [12], the dashed curve (blue) is from the accelerator decay search by 
PS191 [H] (labeled by iV — )■ e"'"e~z/), and the dot-dot-dashed (orange) curve is from the 
peak search in kaon decay [13] (labeled by — )■ fi'^N). 

4 Conclusions 

We have studied the production and the detection of the atmospheric sterile neutrinos. 
With the mass range 1 MeV < M^r ^ 105 MeV, the sterile neutrinos are produced by 
muon or pion decays. The main decay mode of the sterile neutrino N is N Su (88%) 
and the subdominant mode is N ^ e~e^v (12%). Interestingly, the detection via the 
subdominant channel is feasible at Super-Kamiokande. To estimate the amount of the 
sterile neutrino flux (p^i we have reconstructed the parent muon and pion integrated 
fluxes from the well-established active neutrino fluxes 0,^ by making a power-low ansatz 
for the shape of the spectrum. We have calculated the sterile neutrino flux from the 
reconstructed parent fluxes and the energy distributions of each decay mode. The phase 
space suppression of the muon decay /i^ — )■ e^i'e{T^e)N is so effective that the estimated 
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sterile neutrino flux is much less than the naive expectation of (/)Ar = |G^p(/)^+p^ for 
> 40 MeV. 

For the detection, we have calculated the number of the decay event — )■ e~e"*"z/ 
at SK. The upper bound for the muon-type mixing is estimated by comparing the 
invariant mass distribution of the pair and the observational data of the 2e-like rings. 
To estimate the bound, the visible energy cut > 30 MeV is taken into account in 
accordance with the SK data. Moreover, we have studied the opening angle of the 
momenta and clarify the impact of the minimal opening angle under which the two fuzzy 
rings cannot be identified separately. The upper bound for higher mass regime M^v ^ 
50 MeV is significantly relaxed compared to the previous estimation due to the phase 
space suppression of the muon decay. In addtion, the inclusion of the kinematical cuts 
also relaxes the bound for lower mass regime and changes the shape of the bound entirely. 

Fig. [9] shows the upper bound in comparison with the other experiments (see for 
example, Ref. [23l[24]). In the mass window of 34 MeV < < 64 MeV, the atmospheric 
bound is stronger than the peak search in the meson decays. However, the bound from 
the accelerator decay searclill is stronger than the atmospheric bound by a factor of two 
at Mjv = 34 MeV and by a factor of 10 at Mjy = 64 MeV. We nevertheless believe that 
the search for the atmospheric sterile neutrino is an issue of great interest. One reason 
is that not a few progresses in particle physics have been made in cooperation between 
artifical and natural source experiments. Since PS 191 is the sole experiment which set 
the bound stronger than the atmospheric one in the window 34 MeV < Mn < 64 MeV, 
the atmospheric sterile neutrino can play a role of a unique follow-up experiment with a 
natural source. 

Moreover, many points are to be ameliorated for better sensitivities of the atmospheric 
sterile neutrino. It is obvious that our naive scheme to estimate the upper bound does not 
reflect the accumulation of statistics and more advanced analysis should be performued 
for precise argument. The 2e-like ring data for SKI, SKI and SKU [25j can be combined 
to squeeze the room remained for the sterile neutrino. Furthermore, the analysis of the 
2e-like events can be customized aiming for the sterile neutrino detection. For instance, 
one can switch off the visible energy cut of 30 MeV and track the direction of the sum of 
the two ring's momenta, the opening angle between two rings and the visible energy. This 
may dramatically change the sensitivity if the background behaves differently from the 



^ We have included the neutral current contribution to iV — > e e+i^ to apply the bound from Ref. |16j . 
See Ref. [13121. 
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signal in the up/down-going ratio, openging angle distribution, visible energy distribution 
and so on. As for the detection facility, Hyper-Kamiokande has been proposed with the 
fiducial volume 25 time larger than SK and hence the number of events increases by 
the same amount, which improve the bound, roughly speaking, by a factor of five. If the 
new customized analysis upgrades the sensitivity for |0^p by a factor of two for instance, 
Hyper-Kamiokande will improve the bound by one order of magnitude and will reach the 
present accelerator bound. With the refinement of the analysis and the upgrade of the 
facility, there will be a good chance to discover the atmospheric sterile neutrino. 
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A Details for the decay processes 

For the sake of completeness, we present here necessary formulas for the decay processes to 
produce the results in Section |2] and El In this Appendix, we flexibly use the symbol 7 and 
X to denote the gamma factor of the parent particles and (twice of) the daughter energy 
in unit of the parent mass. The electron mass is neglected in the following formulas. 

A.l —>■ e^Vf^N 

This decay process is conducted by the charged- current. In the laboratory frame, the 
decay width F' is 




(18) 



where 



1 - 8r2 + 8r^ - - 24r^ln(r) 



Mm 




(19) 



— 



r = 



7 



m 
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Here En is the muon energy in the laboratory frame. The energy distribution of the sterile 
neutrino N is given by (with x — 2£'jv/m^) 



5'iow(7, x) (1 < 7 < 7cr, X' < X < x") 

gUghi^, X) (1 < 7 < 7cr, X" <X < X+) 

5'high(7, x) (7cr < 7, X- <x <x-^) 

(all others) 



(20) 



where 



x-^2r, = 7(1 + r^) - 7^(1 - r^) , x+ = 7(1 + r^) + 7^(1 - r^) , 



7c 



2r 



(21) 



The functions g^iow and g^high are given by 
2Vx'^ - 4r2 



5'iow(7,3:) 
5high(7,2^) 



3/0 
1 



6/oV7'-l 



(2 - 87^)2;^ + 9(1 + r^)-fx + 4r^(27^ - 5) 
(1 + r^){5 + (367^ - 5Qy + 5/) + 24r27 (3 - 27^) x 

(22) 



-9 (272 - 1) (1 + r^)a;2 + 47 (47^ - 3) J + - yiow(7, x). 



A.2 7r+ ^ //+iV 

The decay width in the laboratory frame is 



Stt 



7 



(23) 



where 



/3/ = V 1 - 2(r^ + r2) + {rl - r^f, tn 



M; 



N 



m,. 



7 



m7r± m7r± m^ri 

The energy distribution of the sterile neutrino is given by (with x ~ 2£Jjv/m^±) 



(24) 



V' dx 



1 1 



(1 < 7, a; < X < x^ 



(all others) 



(25) 



where 



(26) 
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A. 3 N ^ e e^Va (a = /U-) ^) 

This process is conducted by the neutral current if Gg contribution is neghgible. In this 
case, the decay width in the laboratory frame is given by 



JV 



' - — sm t^iy + 2 sm 6w ] - • (27) 



1927r3 V4 ; V7 

In terms of the invariant mass Mfg, the visible energy E^e and the opening angle of e 
momenta Oee-, the differential decay width is written as 



ee? '^v^ COS 

^ee), (28) 

^ N J Jo JxZ J -I 



where 

Mi Em 



and 



^ee = J^, ^^ = JJ^^ ^^M^' =7{^ + Zee)±Vl^-^{l- Zee), (29) 



J(Zee,a;^,COS^ee) = , „ . . T o /1 N^ ' ^^^^ 

a/X^ - 42;eeVX^(l - COsOeeF " 8^ee(l " COSfc'ee)'' 

3z r 1 

K('y,Zee,Xv,COs9ee) = — , — /o + /l COS 6'ee • (31) 

Vj^{l-COSeee){xl-4Zeey\-'' 



Here /o and /i are 



/O = AZee ( 3 + 2(107" - ^)Zee " 4) " 407Zee(l + ^ee)x„ 

+42;ee ((3 - 27^) + 2Zee) Xl + 47(1 + Zee)xl - (1 + Zee)xt, (32) 

/l = 4:Zee ( -1 + 2(27" - l)Zee + ^Z^^) - S^fZeeil + Zee)x^ 

+4 (1 + (272 + \)Zee - Zl) xl - 47(1 + Zee)xl + (1 + Zee)xt. (33) 



By integrating over cos^ee, one finds for instance 

:(1 - Zee){l + 2Zee) (1 < 7> < Zee < 1, < X^, < X+), 



1 dTN _ 1 

Ftv dZeedXy ^ ^"i — \ 

= 2(1 - ^ee)'(l + 2^ee) (0<Zee<l). (34) 



Ftv dZee 

It is now straightforward to obtain ( JT6l) from (IMl) . The analytical expressions for ^, 
distribution are rather lengthy and we do not present them here. 
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